Effects of colored noise on Landau-Zener Transitions: Two and Three-Level Systems 
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We investigate the Landau-Zener transition in two- and three- level systems subject to a classical 
Gaussian noise. Two complementary limits of the noise being fast and slow compared to character- 
istic Landau-Zener tunnel times are discussed. The analytical solution of a density matrix (Bloch) 
equation is given for a long time asymptotic of transition probability. It is demonstrated that the 
transition probability induced/assisted by the fast noise can be obtained through a procedure of 
Bloch 's equation averaging with further reduction it to a master equation. In contrast to the case of 
fast noise, the transition probability for LZ transition induced/assisted by the slow classical noise 
can be obtained by averaging a solution of Bloch's equation over the noise realization. As a result, 
the transition probability is described by the activation Arrhenius law. The approximate solution 
of the Bloch's equation at finite times is written in terms of Fresnel's integrals and interpreted in 
terms of interference pattern. We discuss consequences of a local isomorphism between SU (2) and 
SO(3) groups and connections between Schrodinger and Bloch descriptions of spin dynamics. Based 
on this isomorphism we establish the relations between S = 1/2 and S — 1 transition probabilities 
influenced by the noise. A possibility to use the slow noise as a probe for tunnel time is discussed. 

PACS numbers: 73.40.Gk, 05.40.Ca, 03.65.-w, 02.50.Ey 



I. INTRODUCTION 

The interest to Landau-Zener (LZ) modeP^ is con- 
stantly increasing over last decades as it establishes 
one of the most important fundamental notions in non- 
stationary quantum mechanics. Original LZ model de- 
scribes the probability of transition between two quan- 
tum isolated states coupled by a linearly sweeping exter- 
nal transverse field of a constant amplitude and a time- 
dependent longitudinal field which passes through reso- 
nance with the transition frequency. Even though the 
resulting LZ formula seems to be quite simple, it proved 
to be applicable in analysing the experimental data on 
charge transfer particle collisions^. The model has been 
employed in various studies related to charge transport in 
nanostructures^K^l, Bose-E inste in condensates^* ^, spin 
tunnelling of nanomagnet d 13 * 1 ^ , quantum quenched^ttlZI. 
Especially, LZ becomes a corner stone of theories dis- 
cussing adiabatic quantum computing^N 2 ^ due to a pos- 
sibility to enhance a read out of qubits via the Zener 
flip tunneling^. Such a mechanism has previously been 
implemented for flux qubits 2 ^! and may serve also for in- 
verting spin population by sweeping the system through 
the resonance (rapid passage) in ultracold molecules^. 

In realistic systems however, spins states remain con- 
stantly coupled to their environment. Among various 
mechanisms of dephasing and decoherence of LZ transi- 
tions between Zeeman-split spin states, the coupling of 
two-level system both with phonon bath and nuclear sub- 
system should be mentioned. If the spin-nuclear coupling 
strength is weak enough and the relaxation of the nuclear 
bath is fast, than the nuclear dynamics effects can be re- 
duced to a fast random fielcP^S. In the opposite situa- 



tion, if the nuclear subsystem is slow enough compared to 
characteristic tunnel time, its influence on LZ transition 
can be acc ounted by an effective model of a slow classical 
nois d 29 * 3 ^ . For the simplest LZ scenario, the nuclear dy- 
namics can be considered in the Born-Oppenheimer ap- 
proximation, thus that the spins are presumed isolated 
and transitions activated by an external magnetic field. 
Besides, the noise associated with both a hyperfine and 
dipole fields plays an important role in description of dy- 
namical response nanomagnets 3 ^. 

Experiments with molecular magnets^ revealed the 
presence of hysteresis phenomena in nanoscale molecular 
magnets amount which Mni2 and Fes- LZ transitions at 
the avoided crossing between the Zeeman-split spin levels 
produced by hyperfine interactions have bee n poin ted out 
as responsible of plateaus on hysteresis loop d 31 * 32 - 1 . Hence, 
a number of proposals have been suggested, clarifying the 
effects of nuclear bath, noises and decoherence effects on 
the transition probability in linearly driven systems in- 
cluding two and multi-states systems^H22l_ Several com- 
pact analytic results have been derived to describe these 
effects, namely the Kayanuma's formula^^ for a str ong 
diagonal noise and Pokrovsky-Sinitsyn formulcP^ 2 ^ for 
the coupling to a fast colored noise with off-diagonal com- 
ponents. 

Spin transport processes in magnetic semiconduc- 
tors designs are unavoidably suffer to hyperfine inter- 
actions treated as noise source frustrating spins dur- 
ing transmissions^^]. The common way to protect in- 
formation during the propagation consists on adiabat- 
ically applying an external controlling magnetic field. 
BctthauserF^I have recently presented an alternative ex- 
perimental method to protect spin propagation in spin- 
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transistors including diabatic LZ tunable transitions. In- 
deed, they utilized as a controlling magnetic fields a com- 
bination of a spatially rotating magnetic field B s and a 
homogeneous field B. Thus, the spin states are subject 
to both a constant magnetic field (£>) and a fluctuating 
(Overhauser) field. A theoretical attempt to attack such 
problem has been introduced in Refs.[57|and|25] for two- 
level spin systems by means of fast noise associated with 
random hyperfine interactions. In contrast to it, recent 
experiments on spin polarization of nuclear subsystem via 
time-dependent gate voltage in double quantum dots^HUl 
have shown that the fluctuations of Overhauser field are 
rather slow changing dramatically the properties of LZ 
transition. The "minimal theoretical model" should how- 
ever take into account all low energy two-electron states 
in a dot consisting of three singlet and one triplet {S = I) 
state. 

In this paper we consider the influence of both fast 
and slow classical noise on two- and three- level sys- 
tems. We calculate transition probabilities for the noise- 
induced and noise-assisted processes by using density ma- 
trix (Bloch) equation. The analytical expression for finite 
time probabilities for two- and three- level systems are in- 
terpreted in terms of Fresnel's interference. In addition 
to two standard definition of a tunnel time for LZ tran- 
sition by means of internal or external clock, we discuss 
a possibility to use a noise as yet another probe for the 
LZ time. 

The paper is organized as follows. The Section [TT] is 
devoted to discussion of basic equations for LZ transi- 
tion S — 1/2 derived through Schrodinger and Bloch 



Here, C(t) = [cf \ cf \ C M 
for amplitude probabilities and 

n(t) = e(t) • s 
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approaches. In Section III we discuss the noise-induced 



and noise-assistant LZ transitions in two-level system. 
The classical noise associated with fluctuations of Over- 
hauser field is considered as a colored noise with Gaussian 
realization. Both cases of one and two- component trans- 
verse noise are discussed. The Section IV and V contain 
the key equations for three-level S = 1 system subject 
to both fast and slow classical noise. In Section IVII we 
discuss the LZ transition times defined through internal 
and external clock in the presence of noise. The details 
of derivation are sketched in Appendices. 



II. BASIC RELATIONS FOR TWO-LEVEL 
SYSTEMS 

A. Schrodinger picture, spin-1/2 

The time evolution of J\f states of a quantum mechan- 
ical system with a coherently driven total spin S can 
be described by a system of J\f coupled differential equa- 
tions for the amplitudes c[ S) (t) , C ( 2 S) (t) , ... , C { J ] (i) of the 

(S) . . r*r\. . .ifn .... 



states Vi 5, (X),V4 S) W 



,^\t) (h = iy. 



(2.1) 



is a column vector 



(2.2) 



is the total Hamiltonian of the system, S the total spin 
vector involving all the three generators of the group 
SU(2). 



The relevant aspect of Eq.(2.2) with our aim lies in its 



description of Zeeman splitting of spin states in a linearly 
sweeping external magnetic field. This aspect intimately 
refers to the traditional LZ problem and the functions in 
Eq.(|2.2|) are explicited as follow: 



& x (t) = 2A, Q v (t) = and O z (t) = 2at. (2.3) 

Here, a > is the constant sweep velocity, A the tun- 
neling coupling matrix element between states that we 
assume here as real and lasts from t — — oo to t — oo. 



For the special case of two- levels, the problem Eq.(2.1 ) 



leads to a system of two coupled equations which are 
decoupled to give: 



dz 



iX 



1 



rf2 r (l/2) 



iX + 



1 1 



C[ 1/2) (z)=0, (2.4) 



C { 2 1/2) (z)=0, (2.5) 
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each known as Weber's equations^" where z = 
\[2ute~ 171 1 ^ and A = A 2 /2a. Solutions to these equa- 
tions are computed with respect to the initial conditions. 



For the choice C^ 2 \— oo) 



1 and C. 



(1/2) 



(— oo) = i.e. 



the pa rticl e was initially prepared in the state ^^^(t), 
one has 1 2 : 

and 



7V" 



D_ iX (-ifit), 



(2.6) 



a 



(1/2) 



(t) = A-e^D-ix-ii-ifrt). 



(2.7) 



Here, D n (z) is the parabolic cylinder (Weber's^) func- 
tion, <p a phase factor and fi = \/2ae~~ l7T / 4 . The pa- 
rameter A is introduced hereafter to distinguish between 
the sudden (A <C 1) and the adiabatic (A 3> 1) limits 



of transitions. 



The normalization factors A + and A- 



in Eqs.(2.6) and (2.7) are respectively defined by their 
modulus, \A+\ = \A-\ = y/Xe~ nX / 4 . 

The probability \C { 2 1/2 \t)\ 2 that the system will be 



found in the state ip 2 
given by 



(1/2) 



(t) at any given time t is therefore 



P L z(i)=Ae- 7rA / 2 |£)_ iA -i(-i/it)P 



(2.8) 



Symmetries of levels allow for directly finding the proba- 
bility to remain in the same state. Some asymptotic and 
exact values of Eq.(2.8) are performed with the aid of 
asymptotic series expansion of Weber's functions^. By 
putting for instance t — > oo we recover 

-2ttA 



PlzM = 1 
known as the celebrated LZ formulaP^. 



(2.9) 
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B. Bloch's picture, spin-1/2 



The general solution for time-dependent LZ probabil- 



ity (2.8) is written in terms of products of Weber's func- 



tions. The asymptotic form of this equation casts nev- 
ertheless a very simple exponential form (2.9). In this 



section we present an approximate solution for LZ finite 
time probability (not necessarily in a long time limit) as 
an exponential of single-parametric real function and dis- 
cuss accuracy of this solution for sudden and adiabatic 
limits. 

The time evolution of the quantum dynamical system 



described by the model (2.2) and coherently driven is 



here governed by the von-Neumann equation for the total 
density matrix p{t), 



. dm 

dt 



H(t)J(t) 



(2.10) 



With the help of Eq.(2.10) we find the population dif- 



ference p(t) — p\\{t) — p22 (t) as being solution of the 
differential equation: 



dp(t) 
dt 



-ie-(t)p 21 (t) + ie+(t)pi2(t), 



(2.11) 



where f>2i{t) = Pi2{t) involving p*{t) — p{t) with, 

:f_ x exp (if^ e z (r')dr')e_(ti)p(ti)*i 



2 exp 



(t/! oo e*(r')dr' 



and <d±(t) = Q x (t) ± iQ v {t). Indic es 1 a nd 2 deno te the 
two- level crossing. Inserting Eq.(2.12) into Eq.(2.11) 



with reference to the context of LZ theory i.e. <d z [t) 
2at, without loss of generality, we find an equation 



dp(t) 
dt 



dti cos 



(2.13) 



which can be included in the fami ly of Volte rra's 
integral-differential equations^ 21 . Eq.( |2.10[ )- |2T3l ) cor- 
respond to Bloch's transformations related to the op- 
tical Bloch's^ equation b = —Qxb. Where b is the 
Bloch's vector set on unit sphere by the condition Tip — 
Trp 2 = 1. The z-component of it being expressed as a 
linear combination of diagonal matrix elements of p(t) as 
b z (t) = p(t) = Pll (t) - ^(f) satisfies Eq.fl2~T3l 



d 
dt 



-4A' 



a{e-i{) pW^dtu (2.14) 



for the conventional LZ problem. Superscript (0) refers 
to LZ problem in the absence of noise 



The integral-differential Eq.(2.14) can be solved iter- 
atively with the condition p( ' (— oo) = 1 that preserves 



the total population at any arbitrary time t. A pertur- 
bation series expansion investigation with respect to the 
parameter A <C 1 (see AppendixjA]) is achieved as, 



(t) = 1 + 2 - 2Tr\F(t) + j ( - 2n\F(t) * 



-G(t,X) 



= -1 



exp - 2n\F(t)) - G(t, A) 
(2.15) 

The function F(t) is defined as 



F(t) = 



2a \ In 2 
7T V + 2. 



(2.16) 



c(^j2a/irt) and s{^/2a/irt) are respectively the cosine 
and sine Fresnel's integrals^. The function F(t) is more 
often encountered in the theory of light diffraction where 
it relates to the intensity of light passed through a semi- 
infinite plane bounded by a sharp straight edge with t 
assuming the lateral distance of the edge from the point 
of observation^. 

The function G(t, A) includes all corrections to the ex- 
ponential solution and is determined so that as t — > oo 
asymptotically one comes back to the conventional LZ 
formula Eq.(2.9). Then, for convenience we write our 



finite time transition probability as follow; 



53. 



P-Lzit) = 1 - exp { - 2ttA[F(*) +\nW(t)]}, (2.17) 



(2.12) with 



mw(t) = -— x In 

in which 



1 - 



G(t, A)cxp|27rAF(i)| , (2.18) 



G(t, A) = Xe-^D^i-iz)? (l - e - 2 ^«) . 

(2.19) 

One can see that in the limit A -C 1 the correction 
function lnVT(t) < F(t) for all times. The two forms 
Eq.(2.17l and Eq.( |2.8[ ) are equivalent with the only dif- 
ference being that E q.(|2.17 ) is the exponential form of 
Eq.(|2.8|). The form (12.171) we obtained will be used for 



analytic derivations of finite times transition probabili- 
ties in the limit of slow noise driven LZ transitions. 

The Fresnel's integrals give rise to Fresnel's oscilla- 
tions (see Fig [l] below) and suggest interferences between 
states around the anticrossing region. Fresnel-type oscil- 

_i_ . ,2 

latory factors e involved in Fresnel's integrals origi- 
nate from the phase 



e z (t')dt' = at 2 , 



(2.20) 



accumulated by the two components of the wave function 
during a linear sweep. 
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FIG. 1. Time evolution of the LZ transition probability in 
the diabatic basis of a two-level system. 

The function Plz (t) could have also been found by nu- 
merically solving the time-dependent Schrodinger's equa- 
tion (|2.1|) . A typical result is depicted by Fig[I] above 



where the time evolution of the transition probability 
Phz{t) discriminates two regions. In the domain t < 
all probabilities are smooth monotonic functions of time. 
First, around the crossing time t = 0, one sees a sharp 
transition of Plz (t) which rises from zero to its maximum 
value and later behaves asymptotically with decaying os- 
cillations around the saturation value Plz (00) ■ This last 
performance characterizes the system for small values of 
A (sudden limit) while for largest A (adiabatic limit) oscil- 
lations are strongly mitigated. The first region identified 
characterizes the jump time Tj ump while the second deter- 



mines the relaxation time 



Troiax^- The jump time r jump 
called hereafter LZ transition time is denoted as tlz- In 
adiabatic limit (slow passage) this time should be some 
functions of the velocity a and the splitting energy A 
(see discussion in Section VI I : 

(2.21) 



t lz = C(a, A). 



The rapid drive of a two-level system produces re- 
peated LZ transitions after passing the resonance, drives 
the system into a coherent superposition of states 
where they may interfere and generates Stiickelberg 
fringes^HSni As a consequence, appearing fringes en- 
code informat ion about the system's evolution and en- 
ergy spectrurrP£iiM>l The system will not then feel the 
gap and tlz should be independent on the energy split- 
ting A: 



t lz = C{a). 



(2.22) 



The slow drive in the opposite extreme limit produces os- 
cillations of very weak amplitudes so that one can assume 
a single transition; Stiickelberg fringes of the former sort 
could be absent on the interferometer. The advantage 
being the possibility to probe spectroscopic information 
about the coherent evolution of the system^. As the 
system feels the gap, 



r LZ = £(a,A). 



(2.23) 



Semi classically, tlz is a time necessary to reach a turning 
point on imaginary time axis of integration contouil^. 

This mechanism of slowing down the sweep rate in or- 
der to collect spectroscopic information about a quantum 
system was recently employed by Berns and co-workers in 
their experiment of Spectroscopy analysis of solid-state 
artificial atorrPIl. In that experiment, they pointed out 
that for Stiickelberg interference to occur, the time in- 
terval between two consecutive LZ transitions should be 
much more smaller than the relevant decoherence times. 

An alternative way to tackle the traditional LZ prob- 
lem consists on transforming the equation (2.141 to a dif- 



ferential equation. It can be achieved by applying a sec- 
ond o rder time derivative to both sides of the equation 
(2.14) and excluding integral term with sin [a(t 2 — tf\. 
Hence, we show that (t) satisfies the third order dif- 
ferential equation, 



dT 3 



P i0) (r) 



T dT z I T 



2A)^<»><r) 



= 0, 



(2.24) 



which can be interpreted as a differential equation for 
the probability. Here, we performed the time scaled 
transformation r = ty/a. While amplitudes found from 
the linear Schrodinger's equation satisfy a second order 
linear differential equation, probabilities from the von- 
Neumann equation rather satisfy a third order linear 
differential equation. Similar equation was written in 
Refs.[54j I62H64] in the same context of LZ theory and 
numerically solved in Ref.[S3] with the aid of a Runge- 
Kutta algorithm. Here, using a correspondence between 
Schrodinger and Bloch approaches we propose an analyt- 
ical and exact solution to this kind of equations. 

The natural initial condition p^(-oo) — 1 was grad- 
ually translated at each step of derivations and we solve 



Eq.(2.24| with the conditions: 



dr 2 



,(0) 



(t) \ t =-oo= -8A, 



dr 



P i0 Hr) |r=- 



>=0. 
(2.25) 



As T-L{t) in the Schrodinger's equation (2.1) realizes 
SU(2) symmetry operations, the amplitudes C^^ 2 \t) 

and C2 1 ^ 2 \t) also realize the same set of operations. 
This is due to the temporal linearity of the Schrodinger's 
equation which preserves symmetry properties. While a 
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Schrodinger's equation describes dynamics of wave func- 
tions (SU(2) spinors), the Bloch's equation deals with 
evolution of probability densities combined into a vector 
on unit sphere (Bloch's vector on 2-sphere). However, a 
local isomorphism between SU(2) and SO (3) establishes 
relations between these two objects. 

In the population difference, the occupation probabil- 
ities are express ed in ter ms o f transition amplitudes for 



and p^(r) = | C^^^'* (t) | 2 and this realizes the isomor- 



.(V 



2) 



half-spin in Eqs.(|2.6|) and (127} i.e.: pf^(r) = \c[ 1/2 \t)\ 



phism we talked about. Thus, the solution of Eq.(2.24) 



reads, 



P ^(t) = -A e -^/ 2 [|£_. iA -iHMoT)| 



--\D_ lX (~ifi T)\ 2 \. (2.26) 



Here, /j,q = p/^fa. We may deduce, from the same 
technique, an int egral relation between Weber's func- 
tions. From Eq.(2.12| one may notice that p^(t) = 
Cf /2) (t)cf /2) *(t) and finds 



D_ iX (-i f xt)D*_ i x-i(-ip>t) = -A/x J exp(^[i 2 -£ 2 ])(jL>_ 4A -iH^i)| 2 - I|D_ iA (-i^i)| 2 )dii. (2.27) 



-i (—ifit) can be 



(1/2) *(t)cf /2) (t). Similar 



Similar relation for D*_ iX {—ifjLt)D-i\ 
derived from foi (t) = C\ 

matching procedure was recently employed in Ref. 
to establish an integral relation between Weber's func- 
tion s not from pf^ '(*) as wc did here but from p^(t) in 
Eq.pjl]). 



One can easily check that the limit t^oo applied to 
(2.26) brings us back automatically to Eq.(2.9). Thus, 



the former represents the population difference at any 
given time r. Eq.(2.24), which addressed with LZ tran- 



sitions is valuable to implement an alternative means to 
manipulate the flux qubit and control population differ- 
ence in a micromaser at given time d 65 " 67 ! It might serve 
for transfer of population between two states at any time 
r. For instance by measuring the LZ transition proba- 
bility between two states, it provides information about 
the strength A of the coupling between states. It could 
also offer great advantages in experiments with atoms 
transfer, having only one parameter for control. 



In the domain r < the projection of Bloch's vector on 
z-axis is positive. The system remains in the state where 
it has been set initially. Passing now through the reso- 
nance, p'°)(t) abruptly changes its concavity becoming 
either greater or less than zero. Just around the anti- 
crossing region, the sharp drop of P^°'(t) shows that |1) 
has started to feed |2) via the LZ mechanism. 



In the domain r > 0, one has p^°'(t) < 0, the two-level 
system experiences decaying oscillations while the popu- 
lation difference saturates to finite value. The oscillations 
correspond to an interference between states |1) and |2). 
This last remark tells us that population difference tends 
to maintain the majority of the system into the excited 
state rather than the ground state. 



III. TRANSVERSE NOISE IN SPIN-1/2 
LANDAU-ZENER THEORY 

We now turn into a situation where LZ transitions are 
noise induced. Basically, the coupling between level po- 
sitions fluctuates due to a transverse noise with Gaussian 
realizations. 

LZ effects in presence of transverse classical noise in- 
cluding inter-level transitions are specified by the proto- 



type Hamiltonian (2.2) considering 



Q x (t) = 2f x (t), e y (t) = 2f y (t) and O z (t) 



2at. 
(3.1) 



These definitions are also valid for the case S = 1 
we study below. The mean-zero stochastic functions 
fi(t)(i = x, y) in Eq.(3.I I are characterized by their first 



and second order moments, 

(/*(*)) - 0, (fi(t)fj(t')) = 7?%- exp(- 7 |f - t'\). (3.2) 

Here, r\ stands for the noise intensity that might be re- 
lated to the absolute temperature via the universal fluc- 
tuation dissipation theorenW^ (see detailed discussion be- 
low). Parameter 7 = l/t no i se defines a time scale asso- 
ciated with the noise. Comparison of t no i se with charac- 
teristic time scales of LZ problem gives us a definition of 
fast and slow noise limits. The dynamics of the system 



is governed by Eq.(2.I3) for the occupation difference: 



dp{t) 
dt 



= -4 



cos 



a(t 2 -tl) /+(*)/_ (ti)p(tOdii 



(3.3) 



where f±(t) = f x (t) ± if y (t). The solution of this equa- 
tion is to be averaged over all possible realizations of the 
two-level system (ensemble average) . The result of this 
averaging is different for the two limits of fast and slow 
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A. Fast noise, spin-1/2 

If noise is fast, the characteristic noise time 
tnoise *C tlz, one can average equation Eq.(3.3) di- 



rectly and decouple the product (f + (t)f-(ti)p(ti)) as 
(f + (t)f-(ti))(p(ti)}. The resulting master equation for 
the average (p(t)) gives the conventional equation for the 
transition probability as the average (p(ti)) doesn't really 
change in the exceedingly short time interval £ = t% — t; 



dt 



(p(t)} = -(l(t)(p(t)). 



(3.4) 



Here, the functional Ct(t) = (l[u>(t)] of the frequency 
u)(t) = 2at is defined through the two-time correla- 
tion function K(\t - h\) = (f + (t)f-(ti)) as Cl(t) = 
fiW(t) +n(-)(t), where 



/oc 
exp 
-OO 



±«3(t)£ n(\z\)d£, 



(3.5) 



are the power spectral densities of noise capturing in- 
formation about environmental effects. For the Gaus- 



sian model we considered, Eq.(3.5) result in a Lorentzian 



The frequency ui(t) is antisymmetric t) = —Q(t), the 
Lorentzian spectral density is symmetric in the Fourier 
space leading thus to fK+)(i) = fi (_) (<)- Eq.^SA^ is thus 
readily solved to give: 



(/»(*)) = (/5(-oo))exp 



n(t')dt' 



(3.6) 



and Cl(t) is readily integrated accounting for the spectral 
density. As a result, we obtained a phase i?(t) accumu- 
lated during an interval of time (— oo,t]: 



0(t) 



4 mdt'=^ 



2 ,2a • 

1 H — arctan( — t) 

7T 7 



and 0(£) = 2i?(i). By putting as usual t — oo one gets 



00 - 4tt 



6>(oo) = 2 / n(fW = —K(0). (3. 

a 



So, the transition probability results from Eqs.(3.6) and 



(3.7) or (3.8) as a combination of the initial condition 



p{— oo) = 1 and the conservation law Trp(i) = 1: 



Pfn = 



1 



1 - e- 6 ' 2 



(3.9) 



Her e, w e choose to write 6 = 6{pd) for later convenience. 
Eq.(3.9) is nothing but the Kayanuma's^ result repro- 
duced by Pokrovsk} * 27 ! 28 ! and demonstrating an equal dis- 
tribution of the system between the ground and excited 
states after passing the crossing time for intense noises 
7/ — > oo. (White noise limit can be obtained from (3.2 ) by 
simultaneously taking r/ — > oo and 7 — > 00 while keeping 




-600 - 400-^00 200 400 600 
(a) 



-600-400-200 200 400 600 
(b) 



FIG. 2. Time evolution of the LZ transition probability in 
the diabatic basis of a two-level system in presence of a fast 
transverse noise for the rapid (a) and slow (b) passages. The 
amplitude of the fast noise is fixed. The noise is characterized 
by a dimensionless parameter v = tyt/ 2 /a and dimensionless 
frequency 70 = y/^/a. For all calculations 70 = 100. 



r; 2 /7 = const. Note, that 9 ~ rj 2 and does not depend 
on 7 in that limit.) 

Here, we defined the dimensionless frequency 70 = 
'y/y/ot and dimensionless parameter v = i^rf/a. If 
noise is directed on either the transverse direction (X- 
noise) or the two-components transverse noise (Ay- 



noise), Eq.(3.9) reduces respectively to 



pi 
-if, 



l-exp(-^(/ e (t)/ B (t)))] (3.10) 



and 



pxy _ 1 r 
/" ~ 2 



/ 2ir r 

1 - exp ( - — (f x (t)f x (t)) + (fy(t)f y (t)) 

\ a 1 



(3.11) 



Hence, in order to sum up noises in A and K-directions 
it just suffices to do that in the argument of the expo- 



(3-7) nential in Eq.(3.9). What happens if noise is colored in 



one direction and white in another? The answer to this 
questi on is p rovided by the argument of the exponential 
in Eq.(3.11). Obviously, white noise will dominate the 



colored one and there will not be a complete transfer of 
population: both states remaining constantly coupled. 



The solution (3.7) generalizing (3.8) to finite times 



coincides at very large times with the results obtained 
in Refs.[2"Tl and [2"51. Similar results were discussed in 
Ref.[5S] for nonlinear drive with telegraph noise. A rele- 
vant note similar to the fast telegraph noise in a two-level 
system for this Gaussian model is that, noise fluctuations 
are averaged out as there are no fluctuations (see Figj2|. 

We specify for further purposes that n 2 /a <C 1 and 
rj 2 /a> 1 correspond respectively hereafter to the sudden 
and adiabatic limits of transitions. In the adiabatic limit 
then, the transition probability depends non-analytically 
on the sweep velocity. Thus far, for t no i se -C tlz there 
is no complete transfer of population; the two states are 
constantly occupied. 
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Spin-1/2, in a constant off- diagonal field and a fast 
transverse random field 

Let us consider a spin coupled to a constant off- 
diagonal and sweeping magnetic field and a transverse 
noise source. Such pattern corresponds for instance to 
spin frustrated by hyperfine field or Overhauser field and 
protected by a constant magnetic field. With reference 
to the recent experimental work in Ref.|35] we present 
an alternative way to protect spin propagation in spin- 
transistor (see Introduction). In this frame, LZ tran- 
sitions are noise assisted and noise fields are no longer 
centered at the origin in the A-direction: 



f X (t)=A+f x (t), fy(t) = f y (t). 



(3.12) 



Noise correlation functions for fi(t) are given by Eq.(3.2 1 



The model (3.1) with Eq.(3.12) can also assume a spin 



weakly interacting with an environment, for example a 
nuclear spin bath. Assuming the spin-bath interaction 
as weak enough as bath relaxation is much faster than 
the inverse interaction energy, we may treat the bath as 
a fast noisy magnetic fielcPH 

As far as noise is no longer centred at origin, we call 
p( SF \t) the average of the total density matrix for the 



non-zero mean problem labeled by Eq.(3.12). Straight 



forward calculations for spin-1/2 suggest a governing 
equation of the form: 



-4A 2 



dh cos[a(t 2 - ti)]p (SF) (ti) 



dh cos[a(i 2 - t?)]/4*)/-(*i)P (SF) (*i)- 

(3.13) 

Linear terms of noise have been dropped since after av- 
eraging and use of fast noise requirements they vanish. 
Hereafter, we adopt the label V to denote transition prob- 
abilities related to the non-centered noise. 



As noise is fast, we can readily average Eq.(3.13 1 as we 



did before and apply the decoupling procedure associated 
with the other arguments of fast noise. As a result, a 
differential equation casts a form: 



dt 



<p (sf) W> 



dt lC os[a(t 2 -4)} (p (SF) (ti)) 
-0(i)<p( sp )(*)}. (3.14) 



We obtained the conventional master equation for the 
transition probability in which noise appears as a per- 
turbing source. Noise essentially modifies the standard 
occupation difference p^ (t) by a decaying random phase 
factor. Mainly, the noise produces dephasing during the 
transfer of population. The phase accumulated during 
the short time interval t± — t ~ I/7 is small enough such 
that: 



n(t')dt' 



n(t')dt'. 



(3.15) 




FIG. 3. Dynamics evolution of the population difference in 
presence of a fast transverse noise, 70 = 100 and A = 0.05. 



Indeed, as the characteristic frequency 7 — > 00 it ap- 
pears that ti « t, justifying the approximation Eq.(3.f 5 ). 



Thus, the solution of Eq.(3.f4) can be formally expressed 
as follows: 



(i)) 



exp 



Cl(t')dt' p i0) (t). (3.16) 



For the solution at t — 00 we derive p' ' (00) = 
2e _27rA — 1 from Eq.(2.26 ) and the finding transition prob- 



ability is obtained as follows: 



3 (SF) 

fn 



(3.17) 



It becomes obvious that by setting 9 — we recover the 
LZ formula Eq.(2.9). Likewise the requirement A = 
leads to the solution for fast noise centered at origin. 
Spin state evolution in spin-transistors might be pro- 
tected during the transport by adiabatically applying a 
homogeneous magnetic field (A > 1). This technique was 
employed in Ref . |42) . 

On Fig |3] the abrupt decay of the population difference 
(p( SF \t)) around the anti-crossing region characterizes a 
rapid transfer of population. However, as the condition 
{p (SF) (t)) > is always fulfilled there is no way to expect 
a complete transfer from one of the diabatic states to an- 
other with the fast noise. Fast noise being characterized 
by short time memory essentially creates a dephasing be- 
tween the states of a two-level system. 



B. Slow noise, spin-1/2 

If now noise is slow {t no ise 3> Tlz), the decoupling pro- 
cedure is not applicable, thus the density matrix equation 



of motion (3.3) cannot be reduced to a master equation. 



Instead one has to solve Eq.(3.3) and perform ensemble 



average over the distribution Q of noise: 



(...) = 



f 



2TT1] 



L dQ - exp (-$) 



(3.18) 
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The brackets (...) indicate as usual the ensemble average. 
In a given realization of classical field Q the LZ proba- 
bility is given by standard equation 



Plz(Q) = 1-exp 



(3.19) 



If the noise in transverse direction is single-component 
we can always rotate our coordinate frame such a way 
that the fluctuations occur along X-direction. The noise 
- averaged LZ probability is defined as 



1 f°° / x 2 \ 

and after straightforward calculation is given by 
P x = 1 - 1 

sn 



1 



27T77 2 



(3.20) 



(3.21) 



If the transverse noise is described by two orthogonal 
non-correlated components, the transition probability is 
averaged with a two-dimensional Gaussian distribution: 



P ™ = 2^rfJ dX j d y p Lz(x,?;)exp( 



which after calculations acquires the form 



x 2 + y 2 



pxy 
sn 



1 



1 



2TT7] 2 



2ri 2 
(3.22) 



(3.23) 
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FIG. 4. Time evolution of the LZ transition probability in the 
diabatic basis of a two-level system in presence of slow one- 
and two- component transverse noises (see discussion in the 
text). 



In the limit t — > oo, where F (oo) = 1 and lnlT(oo) — > 
we return to Eq.(3.21|. One can do the same with the 



two-components transverse noise (XY"-noise model) and 
get 



The difference between Eqs.(3.21) and Eq.(3.23l is a 
consequence of the effective two-dimensional character of 
noise fluctuation spectrum in the latter case and its one- 
dimensionality in the former case (see alscE^). 

As expected, renormalization of the inter-level distance 
by a stochastic function considerably affects the generic 
picture of LZ transitions with off-diagonal coupling. Fast 
noise pours a large energy into the system during the 
crossing, destroys the memory of the system in the do- 
main of strong couplings identically distributing the sys- 
tem between the ground and excited states. We note, 
that in contrast to effects of the fast noise which does 
not change analytical properties of the LZ probability at 
extreme adiabatic limit a — > 0, the two-component slow 
noise transfers the property of the LZ probability from 
Gaussian to Lorentzian, thus making it analytic function 
of a in this limit. 

The expansion method we exposed may allow one to 
formulate transition probabilities for finite times. For the 
case of X-noise for example, solution in Eq.(2.17l gives: 



1 



(3.24) 



F(t)+lnW(t) 



pfm = i - 



2-7TT7 2 



1 



F(t) + \nW(t) 



(3.25) 



Similarly the limit t — > oo brings us back to Eq.( |3.23| . 
These solutions may now be interpreted in an interference 
pattern as they involve Fresnel-like integrals via F(t). 
Pg n and P*y have the same shape but Pgv is greater 
than P* ■ 



Spin-1/2, in a constant off- diagonal field and a slow 
transverse random field 

We dealt up to point with LZ transitions induced 
by the slow noise. This approach can be generalized 
straightforwardly to the case where LZ transitions are 
induced by an external magnetic field (in the case of spin 
systems) or by an effective "field" associated with the 
finite transparency of the inter- well barrier in a double- 
well potential in cold gases. To describe this effect, one 
should take into account the value of the stochastic field 



fi(t) as defined in Eq.(3.12|. Then, averaging over slow 
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one-component noise fluctuation, results in 
cxp(-27rA$i(£)) 



-Psn(t) = 1 



1 



2ttt] z 



(3.26) 



F{t) + \nW(t) 



where the phase $> v (t) with v = 1,2, expressed as 



1 



+lnW(£) 

27TW77 2 



+ In W(t) 



(3.27) 



is due the local deviation of noise created in the X- 
direction. If noise was also shifted along the Indirection, 
we would have an additional phase such that the argu- 



ment of the exponential in Eq.(3.26) would be 3> v (t) = 
&v(t) + ^«(*)- This would offer an access to sum noises. 
Nonetheless, the choice we adopted has a great advantage 
in controlling noise fluctuations in a two-level system. 
For the meantime, the infinite time limit of Eq.(3.26| 
suggests that: 



^n(oo) = 1 - 



1 



27T?f 



: exp • 



2ir\ 
1 + 2nr) 2 /a 



In the sudden limit, 

V^(oo) « PlzM 
In the adiabatic limit, 



r^ n (oo) = i - 



27T77" 



exp 



I 2n 2 l 



}• 

(3.28) 
(3.29) 

(3.30) 



Thus, the argument of the exponent does not depend on 
the velocity. Similarly, slow XY-noise in the presence of 
a constant magnetic field results in 



v:i{t) = i - 

Therefore, 

«°o) = 1 - 



cxp(-27rA$i(i)) 



2tti] 



2 



O 



F(t) + \nW(t) 



(3.31) 



1 + 2irr] 2 /a 



exp 



f 2tt\ -i 

I ~ l + 27ra 2 /a/' 



2irrj 2 1 a. . 



and in the sudden limit, 
while the adiabatic limit reads: 



^(oo) = 1 



I 2t7 2 I 



2-KT\- 



exp 



2?7 2 



(3.32) 
(3.33) 

(3.34) 



A two-level system undergoing simultaneously the effects 
of a slow noise in the presence of magnetic field in the 



sudden limit (rj 2 /a <C 1) seems not be sensitive to noise 
structure for long time asymptotic of transition proba- 
bility. For such a setup, the magnetic field effects prevail 
on the noise and a pre-exponential factor is close to one. 
These effects are supported by Figj5] where we plotted 
Eq.(3.28l for t = oo. We essentially show on Fig[5]that 
adiabatic addition of noise, considerably suppressed the 
previous tendency. 

Besides, by putting rj — > in Eq.(3.28) one imme- 



diately comes back to (2.171, the LZ formula for finite 
times. 

The averaging procedure described in details allows to 
calculate the LZ probability if, for example, the noise is 
fast in one of transverse directions and slow in another 
orthogonal direction. In that case one should first aver- 
age the Bloch's equation over fast realization and after 
solve the "effective" Bloch's equation in given realization 
of slow fields. As a result, the fast noise contributes only 
to the argument of LZ exponent, while the slow noise 
appears both in exponential and pre-exponential factors. 
Therefore, numerical fit of experimental datcP^^ could 
provide an information for both kind of noises without 
requiring additional measurements. 

By applying the fluctuation dissipation theorem one 
can associate some effective temperature with equal time 
two-point correlation function as follows 



(Mt) 2 ) =A-T, 



(3.35) 



where A and T are respectively the coupling constant 
with the environment (model dependent) and the abso- 
lute tem peratu re in the units fcs = 1. In this frame, the 
solution (3.34) acquires the Arrhenius^l form 



a 



2-kA-T 



exp 



{-f} 



(3.36) 



where E — A 2 /2A is the activation energy^!. The pre- 
cxponcntial factor provides proper normalization of the 
distribution. 
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FIG. 5. (Color online) Landau-Zener transitions probabilities 
for two-level system at infinite time as a function of dimen- 
sionless parameter irA 2 /a in the presence of noise for different 
values of v — irrj 2 /a. 
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The theory of noise induced LZ effect may be extended 
to the multilevel LZ problems, where more complicated 
patterns of transient oscillations in the tunneling proba- 
bility of transition from initial to final state of a nanosys- 
tem with nontrivial dynamical symmetry is expected. 



IV. BASIC RELATIONS FOR THREE-LEVEL 
SYSTEMS 

A. Schrodinger picture, spin-1 



The Eq.(4.3) can also be represented in the integral- 

1 



differential form: 



dz 



C^\z) = -4iA 



dz\ cosh 



-(z*-z() C£\ Z1 ) 



(4.4) 



reminiscent to (2.14 1. In order to find the solutions, it 



would be more instructive to realize that the square of 



Weber's functions in Eqs.(2.6) and (2.7) satisfy the third 



order differential equations considered by Mitra in Rcf. 



Let us consider the LZ transition in a three-level 
(spin S = 1) system. The upper and lower levels 

denoted as cj 1 ^) and C^fy) are characterized by 
the eigenenergies E±{t) = ±2(a 2 t 2 + A 2 ) 1 / 2 as in the 
LZ process. A middle level wavefunction is denoted 
by C^it). Corresponding eigenenergy E(t) = does 
not evolve in time while the transitions between only 
neighbouring levels is allowed. 



Let us define the diabatic states. The projections — 1, 
and +1 of the spin are to be matched to the states |1), 
|0) and 1 2) respectively to form the only avoided-level 
crossing point. 



In the background of equation (2.1) and the given 
parameters in (2.3), one obtains a system of three 



decoupled differentials equations for the states c[ (t), 

C$~\t) and C*2 (i). The first pair of equations for the 
states with minimal/maximal projection to z-axis is 
given by: 



^f»(z) + (4zA - 2 - *)±C?\z) - *cf >(*) = 0, 



(4.1) 



^C^(z) + (4zA + 2 - *)±C?\z) zC^(z) = 0. 

(4.2) 

The dynamics of the middle level is independently de- 
rived and is governed by a third order linear differential 



equation of the form (2.24). For sake of consistency it is 



presented here as follows 



= 0. 
(4.3) 



When obtaining these equations, no specific assumptions 
on preparation of the system have been adopted. One 
may arbitrary select a particular ground state and ask 
questions about the probability to find the system on 
the excited state. 



dz 3 



C[ 1/2 \z) 



(4zA 



z C{ 1/2) (z) =0, (4.5) 



dz 3 



a 



(1/2), 



(z) + (4zA + 2-z 2 ) 



dz 



a 



(1/2: 



C, 



(1/2), 



(z) =0. (4.6) 



The solutions to these equations are, D 2 x (z), D 2 _ iX _ 1 (iz), 



D 



_ iA _ 1 (— iz) or any products of two of D 2 x (i 



D_ iX Aiz), and D_ iX _ 1 (~iz). The solution of Eqs.(4.1 ) 



and (4.2) with given initial conditions can therefore be 



constructed with the basis of these functions. 

We first consider the initial conditions C[ (— 00) = 1, 



Cq(-oo) = and C^'(—oo) = 0. For these choice, 



the solution read, c\^\t) 



B_D 2 _ lX (-zz), Cf(t) 



B + D 2 _ iX _ 1 (—iz). Here, B± are normalization factors 
B± = \A±\ 2 . Basically, we establish relations between 
the wave functions of triplet S — 1 state and spinor 
5 = 1/2 states: 



C[ 1) (t)= C? m {t) 



.(1/2), 



a 



(1/2) 



(*) 



(4.7) 



(4.8) 



(4.9) 



We solved (4.3) with the aid of the normal izat ion condi- 
tion 8 T,m=-s\ C m } (t)\ 2 = 1 and obtained JZT 



This matching procedure could be applied to the sec- 



ond special case C^(— 00) = , Cq X> (—oo) — 1, and 



(i)/ 



C^j—oo) — since with the condition Cq~'(— 00) = 1 



(i), 



Eq.(4.3) or its integral-differential representation (4.4) is 



structurally isomorphic to (2.24) and 



C£'(t) = -\e-* x l 2 \D 



-iX-\{-i(it)\ 
1 



A 



|D_ iA (-^)l • (4-10) 
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In order to satisfy the conditions c[ X \— oo) = 
C^\— oo) = one should construct solutions which con- 
tain D-i\-i(— iz). One can explicitly check that the so- 
lutions of Eqs.(4.1) and (4.2 1 satisfy the identities: 



C ( 2 \t) = V2Xe-^ /4 



1 



' , M>l k 2[* 2 - a ?]J C 1) (*i) <fe ii 
(4.11) 



exp 



(4.12) 



and 



d 

dz 



C ( 1] (z) = >/2A(cf + ^ 1} (z))e-"/ 4 . (4.13) 



Relations (4.11) and 



representations of (4.1 



4.12| turn out to be the integral 



and (4.2) 



Inserting the exact expression of Cq (t) given in 
Eq.(|4.10| into Eqs.(|4.11|) and (|4.12|) and using identity 



(2.27) one obtains after straightforward calculations: 



cf>{t) = \/2Aexp(V 



ITT 7rA s 

D^i-ifxt^D^ti-ifit)]*, (4.14) 



2A exp 



(iip' 



in 7rA N 

[U_a(-t/<t)]*I>-u-i(-*/*t). (4.15) 



Here also represents an arbitrary phase factor. An ex- 
tension of the conventional LZ theory 5 = 1/2 is there- 
fore achieved for this case via the connecting relations: 



c£\t) = \c^Xt)\*-\c^\t)\\ 



(4.16) 



(4.17) 



(4.18) 



For completeness, the solution of (4.1 ) - ( |4.2[ ) with initial 
conditions (Cf } (-oo) = 0, C^j-oo) = 0, Cf \-oo) 



1) can be found with the help of (4.16 )-(refequa54). We 
would like to point it out that the three-level 5=1 
problem possesses an additional level symmetry corre- 
sponding to a particle-hole symmetry of SU (2) group^. 
A time evolution operator for 5 = 1 is then constructed 
as follows: 



cf /2) (i) 



V2c[ 1/2 \t)d 2 1/2) (t) 

-V2C[ 1/2 \t)4 1/2) *(t) |Cf /2) (i)| 2 - \C { 2 1/2) (t)f 



a 



(1/2), 



(*) 



-V2C[ 1/2 >(t)4 1/2 >(t) 



a 



(1/2) 



(*) 



/2C[ 1/2) *(t)C^ 1/2 \t) 



C[ 1/2) *(t) 



(4.19) 



The time evolution operator in (4.19) is thus a general- 
ization of that obtained in Ref . [2S\. For comparison with 
[28j we present all transition matrix elements in the limit 
t = oo in the TABLE|U 

We have numerically computed the dynamical evo- 
lution of the model (2.2) for 5 = 1 by solving the 



Schrodinger's equation for the amplitudes Cn\t), (n = 
1, 0, 2) and plotted the populations \Cn\t)\ 2 for the sec- 
ond projection of 5 we adopted above. The results are 
depicted on Figj^j As expected, populations |C^(t)| 2 
and | C^ 1 '* (t) | 2 are identically distributed on levels |1) and 
|2> • 

These figures also show how population of each of the 
levels changes with the LZ parameter. Analytically, one 
recovers this dynamics by taking the limit t — > oo. The 
results perfectly agree with the CarolP^l results (see TA- 



BLE j!]). 

We focus hereafter exclusively on Poi(t) = \^oi{t)\ 2 
and -Poo(i) = l^oo (i)| 2 deduced from the second block 
of the table which corresponds to the second preparation 
we discussed above. 



B. Bloch's picture, spin-1 

The dynamics of the system is here governed by 9 equa- 
tions for the components of the 3x3 density matrix which 
might be reduced to a set of 6 equations due to symmetry 
of the levels. The condition for the probability conserva- 
tion tells us that the knowledge of two of these elements 
is enough to compute the remaining. 

In this regard, we reduce the problem to a set of two 
equations for the differences p+(t) = pu(t) — poo(t) and 
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Initial Occupation t 



-oo | Final Occupation t = oo 



2-n-A _ g — 47rAj 
2ttAi2 



2[e -"" - f 
fl-c 



2[e" 

j 1 _ 2e -2^A]2 

2[e" 



fl-e" 



2[e" 



TABLE I. Landau-Zener transition probabilities in a three- 
level system. 



that the Schrodinger's dynamics of TV-level system de- 
scribing spin 5 = (TV — l)/2 is given by set of TV coupled 




(a) 



(b) 




FIG. 6. Time evolution of the LZ transition probability in 
the diabatic basis of a three-level system. 



P—(t) — f>22{t) — Poo(t) where the indices 1, and 2 enu- 
merate the three crossing levels. 

It would be convenient to express the probabilities 
Poi(t) and Poo(t) in exponential form as we did in the 
preceding section for LZ spin-1/2 probabilities. From 
Eq.(4.10) a simple connection between the occupation 



probability Poo(i) and Plz(^) is deduced as Poo(i) = 
(2Plz(£) — I) 2 - Then, the normalization of probabilities 
help one to express the finite time tunneling probabilities 



Poi(t) = 2 cxpj - 27rA[P(f) + lnW(*)]}- 
exp j - An\[F(t) +lnW(t)}\ 



(4.20) 



and 



Poo(*)= 1 - 2 cxpj - 2wX[F(t) + lnIU(t)]} 2 , (4.21) 

which will help us later to derive the LZ transition prob- 
abilities in the slow noise limit. 

The conventional notation adopted throughout this pa- 
per will be reflected by the probabilities restricted to the 
two limits of fast and slow noise. Detailed calculations 
for relevant equations are deferred in Appendix [Bj 

To summarize the discussion on correspondence be- 
tween Schrodinger and Bloch approaches we point it out 



first order linear differential equations (LDU) resulting in 
TV independent TV-th order LDU. On the other hand, the 
time evolution operator could be constructed with the ba- 
sis of fundamental spinor representation of SU(2) group 
by means of Jacoby TV — 1 order polynomials (see also 
Refs. frjTlITT] '). Therefore, in general the solution of TV-th 
order LDU is expressed in terms of TV — 1 power of We- 
ber's function. The Bloch's dynamics of spin 5, however 
is based on one vector and 25—1 tensor Bloch's equa- 
tions since the density matrix has 25 conserving values^. 
The connections between solution of integral-differential 
Bloch's equations and LDU Schrodinger's equations es- 
tablish integral relations between products of Weber's 
functions. 



V. TRANSVERSE NOISE IN SPIN-1 
LANDAU-ZENER THEORY 

In this section we evaluate the tunneling probabilities 
for the case when the inter-level distance between the 
states of a three-level system is renormalized by a ran- 
dom classical field. Similar studies were performed by 
PokrovskjPS with restriction to fast noise. To our best 
knowledge the slow passages have not yet been investi- 
gated. For this purpose, we still use the argument ac- 
cording to which as noise is slow, the LZ transition prob- 
abilities can be obtained from an ensemble averaging over 
all possible noise realizations. 

The equations (Bl )-(B5) combined together in a fash- 
ionable way read: 
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dt 



= -4 



e + (£)e_(ii)cos a{t 2 -t\) p+(h)dt 



6 + (f)e_(ti)cos a(t 2 -tf) p-{ti)dtx, (5.1) 



dp-it) 
dt 



6 + (t)e_(ti)cos a{t 2 -t\) p+(ti)dt 



9 + (t)6_(ti)cos a{t 2 -t 2 ) p-(h)dt 



(5.2) 



Therein, 9±(i) = (6 x (i) ± iO v (t))/V2. As the func- 
tions 0+(t) and O-(t) involved noise fields, we ignored 
all terms containing products of the forms 0+(t)0+(ti) 
or e_(t)9_(ti) and obtained Eqs.fO) and $3. The 



reason is simple; in principle, all averages (0 + (Fj8+(£i)) 
or (0_(i)0_(ii)) do not contribute as well as they vanish 
as Gaussian correlators. 

Obviously, this procedure will not be applied to the 
case where 0+(t) and O-(t) are not noisily correlated. 



As can be seen from Appendix [B| the components p\2{t) 
or P2i{t) would be involved. 

We now check whether the method we employed in Sec- 
tion III for conventional spin-1/2 LZ transition might be 
related to a third order differential equation. As we just 
explained in the preceding paragraph, additional terms 
relat ed to p\2{t) or p2i{t) should enter Eqs.(5.1) and 
(5.2). One has as a result: 



dt 



pf{t) 



dt 



-8A 2 



a(f-ii) p { +\ti)dtx -4A 2 



-6A 2 



p ( -°\t) = -4A 2 











f 


exp 


ia(t 2 


~t\) 


> — oo 








COS 


a(t 2 


~t\) 





a(t 2 -t'i) p { °\h)dt 



exp 



ia(i l -i{) p^(tx)dtx, 



8A^ 



COS 



-6A Z 



exp 



ia(t 2 -tf) p { £>(tx)dtx -6A 



a(t 2 - t\) 



P { -\tx)dh 



exp 



ia(t 2 - tf) p^(ti)dt 



(5.3) 



(5.4) 



It is less obvious to derive a third order differential equa- 
tion for the function (p?(t),p_(t)) as we did for spin 
S = 1/2. The functio n [p 2x{t), Pi2( t)) generated by the 
last two terms in Eqs.(5.3) and (5.4) might be viewed as 



external sources for an homogeneous matrix equation of 
the form Eq.(2.14). However, an appropriate choice of 



variables leads to the non homogeneous equation: 
d 



dt 



5<°)(t) = -4A 



cos 



a{f-ti) p {0) {tx)dt 



-6A^(t)Q(°)(t). (5.5) 



Obtaining this, we define corresponding functions in the 
absence of noise as follows: 



5<°>(t) 



pf{t) 
P { °\t) 



and Q (0) (t) 



P^{t) 



(5.6) 



and the square inter-level distance is given by the matrix 



form 



A^=2A 2 



1 1/2 X 
1/2 1 , 



(5.7) 



The operator J{t) acts onto the subspace Sb of t he v ector 
Q^ (t) and generates the last two terms in Eqs.(5.3) and 
( |5.4[ ) . The Lie algebra associated with the time derivative 
properties of J{t) can be defined from the relation, 



J(t)p { 1 °2 ) (t) = J COS 

sin 



a(t 2 -tj) Rc[p[ 2 ] {t^dt! 
a(t 2 - tf) Im^Oi)]*!' (5 



Eq.(5.8) shows J{t) in action into Sb- The action of 

J(t) is symmetric so that J{t)p${t) = J{t)p${t). Ob- 
viously, the operators d/dt and its higher orders d q /dt q 
(q > 1) have non-zero actions in Sb ■ The simultaneous 
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actions 4 o jf(i) and jj^ o jf(i) including 4^ o J~(i) are 
then also defined in 1S5. 

By evaluating f t J{t)p { ^{t) and §zJ(t)p${t) and 

considering the initial condition P12 ( — °°) = we define 
equivalent initial conditions helpful for further purposes. 
It can be verified that, 



J(t)p^(t) | t =-co= J(t)p^>(t) \ t =-oo= 0. (5.9) 



dt 



(0), 



oo) = —(1 1) T . The solution is given by 



P+ 



(r) = Ae-' rA (2|Z?_. iA (-^or)| 2 |C- l A-i(-^or)| : 



+A 



|-D- 4 a-i(-«A*ot) 



(5.13) 



The nontrivial dynamics of population difference for 
three-level LZ transition at any given time r is governed 
by Eq.(5.10|. Considering the limit r — ¥ oo, one finds for 



the population difference 



These properties are directly applicable to pt$ ' (t) as 
P21 (~ 00 ) = 0- We omitted the symbol o keeping in 
mind the ordering of actions in S^, namely, j{t) first 
passes and d q /dt q follows. 



The Eq.(|5_5|) mimics Eq.( |2.14[ ) for the vector p w {t) 
function of one variable. It may be transformed to a non 
homogeneous linear third order differential equation for 
a vector function: 



d?_ 
dr 3 



P {0) (r)- 



1 d 2 



t dr 2 
- 2 ' 2A 



^°)(r)-4 
d 



2A 



M . 



P i0) (r) 



M) Tr P 



(r) =-6A 2 X(t), 



(5.10) 



pf (00) = p w (00) = 12e- 37rA sinh(TrA) - 1. (5.14) 



Using Eq.(5.14) and the conservation of probability 



Trp(r) = 1, one recovers the transition probabilities in 
the second block of TABLED 



A. Fast noise, spin-1 

For the convenience of reader we review briefly the dis- 
cussion about the effects of fast noise on a three-level sys- 
tem. In the spirit of previous derivations, we transform 



Eqs.fjSJjJ-lJS^ to 

d 



dt 



Pit) 



-4 I cos 

-00 



a(t 2 - tl)\n M (\t - hDpit^dh. 

(5.15) 



where 



X(t) 



J a dr z a dr 



4r' 



Here the vector 

P(t) = 

and the matrix correlator 

n M (\t-t!\) = 2 



(p + (t)) 
(p-(t)) 



1 1/2 
1/2 1 



(5.16) 



(5.17) 



with Am = A M /2a standing obviously for a matrix LZ 
parameter. The best way to see X(t) is to consider 
pfkr) = C?\t)C^{t) and p^(r) = cf * [t)C? ( T ) 
which in their primitive sense would include Weber's 
functions. To solve these equations we fix the initial con- 
ditions for the system set in the state ip^ (r) as done in 
Section IV A Accounting for Eq.(5.9), we solve Eq.(5.10) 
with 



^M ] (t) as 



are defined^ The equation (5.15) is structurally identi 



cal to Eq.(3.3). Similarly we define Cl M (t) = Cl$(t) + 



^(*) = 



exp 



±iQ(t)SK M (\€\)d£. (5.18) 



The solution of Eq.(5.15) can be found with the help of 



the matrix transformation 



dT 2 ? 



(<>) 



(r) 



-8A 



M) 



The adjacent condition f>^(— 00) is fully defined as and reads as follows: 



dr 



P {0) (r) | T =-oo=0. 



(5.12) 



exp 



n M (t')dt' 



cosh - — sinh 

6 6 

— sinh - cosh - 

2 2 



(5.19) 
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(Pn(oo)) = -(l + {p+(-oo)) 
(poo(oo)) = ^(l - (p+(-oo)) 



sinh — h 2 cosh - 
2 2 



cosh sinh - 

2 2 



(p 22 (oo)) = - (1 - (p + (-oo)} 



cosh 



2 sinh 



Initial Occupation t = — oo | Final Occupation t — oo 



3 p -S/2 
2 



1(1 



-36>/2\ 



-36/2\ 



1(1 



1(1- 
|(1 + 2e- 39 / 2 ) 



1(1 



-36>/2\ 



1(1 



3 -9/2 , 1 - 
2 ~ 2 



3S/2\ 



1(1 



-39/2\ 



l(l + | e -8/2+ l e -39/2 



TABLE II. Fast noise transition probabilities in a three-level 
system. 



These results agree with those of Pokrovsky^S The 
general form of these equations for arbitrary t can be 
obtained by replacement 8 — > 9(t) (see TABLE|ll| where 
6{t) defined similarly to On TABLE|IT) we show 

infinite time transition probabilities for all possible initial 
preparations of the system. 

One can see that the transition probabilities for S = 1 
have the same structure as those for S = 1/2, namely, 
in the white noise limit one recovers equal probability 
distribution. 



e - e -(/M-co)) 
+ (P-(-oo)) 
e" + <£_(-«>)) 



cosh — h 2 sinh - 
2 2. 



sinh cosh - 

2 2 



sinh — h 2 cosh - 
2 2 



(5.20) 
(5.21) 
(5.22) 



r 



Spvn-1 in a constant off-diagonal field and a fast transverse 
random field 

Let us now investigate the case of the fast-noise as- 
sistant LZ transition. The two-component noise is de- 
fined by Eq.(3.12). Since the mean- value of stochastic 
function describing noise in A"-direction is non-zero, the 
(6±(*)9±(*i)) ^ 0, and the sources induced by p\\{t) 
and pi2 (t) enter the equation for p+ F ^ (i) and p^ F ^ (t) : 



(pi SF) W> 



and Q (SF) (t) = 



(5.23) 



The dynamics of the system is described by equation: 



dt 



3(SF)( t) = _ 4A 



cos 



S(SF) 



{h)dh 



cos 



a{t 2 -t\) 

a^-tl^KMQt-hDp^ih) 

~6A 2 J{t)Q {SF) (t)dt 1 . 

(5.24) 



The requirements with the help of which we tackled prob- 
lems of fast noise so far, proves to be well applicable and 
the dephasing enters the dynamics as the choices 



p( SF )(i) = exp 



and 



USF) 



(t) = exp 



h M (t')dt' /5 (0) (t), (5.25) 



Cl M (t')dt' Q (0) (t), (5.26) 



prove to entirely solve Eq.(5.24). We present solutions in 
general form: 



(pf 1 F) (oo)> = ^(l+pi° ) (oo) 



sinh — h 2 cosh - 
2 2 



(4 S F) (oo)) = ^l-pV U; (oo) 



5(0), 



cosh sinh - 

2 2 



e- e -p W (oo) 



P i0) (oo) 



cosh — h 2 sinh - 
2 2 



sinh cosh - 

2 2 



(5.27) 
(5.28) 
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Initial Occupation t — — oo | Final Occupation t — oo 

1 |(i _ [| e -»/a _ _ 4e -2^ + 3e -4.A ] + [| e -«/a + i e -3fl/2p c -4*rA _ 2e -a*A]) 

1(1 _ e - 3S / 2 [l - 4e- 27rA + 3e- 4 " x ] - e" 39 / 2 ^" 4 ^ - 2 e - 27rA ]) 

|(1 + [f e-"/ 2 + |e- 39 / 2 ][l - 4e~ 2 - A + 3 e - 4 ^ A ] - [fe^ 2 - f e - 3e / 2 ][3 e - 4 ^ - 2 e - 2 ^ A ]) 

|(l + e- 3e/2 [6e- 2,rA -6e- 4nX - 1]) 

1 i(l_ 2e - M / 2 [6e- 2 ' rA -6e- 4 ' rA -l]) 
i(l + e - 3 "/ 2 [6e- 2 ' rA -6 e - 4 ' rA -l]) 

|(l + [f e - 9 / 2 + | e - 39 / 2 ][l - 4e- 2 " A + 3e- 4 " A ] - [fe" 8 / 2 - fe" 36 / 2 ] [3 e - 4 " A - 2 e - 2 " A ]) 

1(1 _ e ~ 3e / 2 [l - 4e- 27rA + 3e- 4,rA ] - e" 39 / 2 ^" 4 ^ - 2 e - 27rA ]) 

1 |(i _ [| e -"/ 2 _ i e - 3 "/ 2 ][l - 4e- 2 " A + 3e- 4 " A ] + [fe^ 2 + § e - 3e / 2 ] [3e- 4 " A - 2 e - 2 ^ A ]) 



TABLE III. Transition probabilities for a three-level system in both a constant magnetic field and a classical transverse noise. 



(pi S 2 F) (oo)> = i(l-pf (oo) 



cosh — + 2 sinh ■ 



£ 0) (oo) 



sinh cosh - 



(5.29) 



The differences (oo) and p_ (oo) can be found from 
TABLE|l] and corresponding transition probabilities are 
presented in the TABLE JTTT| Let's label the three tables 
presented so far respectively by Tl, T2 and T3. Now, 
assume that T3 is a function of and A, i.e.: T3 = 
T3(0, A). It can then be seen that T3(0, A) = T1(A) and 
T3(0,O) = T2(0). 



transverse noise are given by 



P* n [0^1](i) = 2 



cxp(-27rA$i(i)) 



F(t) +]nW(t) 



cxp(-47rA<I>2(t)) 



1 + 



2 r 



F(t)+]nW(t) 

(5.30) 



Slow noise, spin-1 



The transition probabilities for S = 1 subject to the 
slow noise are obtained in the same spirit as was discussed 
for S = 1/2. In a given realization Q of noise, the sys- 



tem ofcquations for the pop ulatio n differences (Eq.(5.3|- 
Eq.(5.4|) is reduced to Eq. (5.10|. The solutions of this 
problem for t = oo are derived via Eq.(5.14). 



P£,[0-X)](t) = l- 



4exp(-27rA$i(t)) 



2nrj 



2 r 



F(t)+\nW(t) 



4exp(-47rA$ 2 (*)) 



F(t) + \nW(t) 

(5.31) 



The LZ solutions for the case of one-component slow 



For the two-component transverse noise, the transition 
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probabilities read: 
W-H](t) = 2 



exp(-27rA$i(i)) 



2irif 



F(t)+lnW(t) 

exp{-4n\$ 2 {t)) 

47T77 2 r 

' F(t)+lnW(t) 



(5.32) 



and 



CT)->0](t) = l 



4exp(-27rA$i(£)) 



2ttt7^ 



- + 



a 



+lnW(*) 
4exp(-47rA<J> 2 (i)) 

F(t)+hxW(t) 

(5.33) 
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FIG. 7. Time evolution of the LZ transition probability in 
the diabatic basis of a three-level system in presence of slow 
one- and two- component transverse noises. 



Note that T e \0 -> 1] = Pf„[0 -> 2], with £ = 1,353/. 



The solutions ( 5.30 )-( 5.33) represent general solution of 
LZ problem for three-level system in the presence of slow 
one- and two- directional transverse noise (both noise- 
induced and noise-assistant). The solution of slow-noisc- 
induced LZ transition can be extracted from (5.30 )-( 5.33 ) 
by putting A = 0. 



VI. DISCUSSION ON EFFECTS OF NOISE ON 
LANDAU ZENER TIMES 

It is well-known that if one deals with a system of con- 
sequent Landau-Zener transitions, it is not sufficient to 
characterize a behavior of such system by an asymptotic 
values of pro babilities. One also needs to define a tun- 
nel tim o 1 54 ! 7 ^ "-^ in order to put a borderline between two 
cases when the consequent tunnel processes can or can 
not be considered independently. 

There exist several ways to define the tunnel Landau- 
Zener time for two level system. Although we are not 
going to dwell onto detailed discussion of tunnel times 
in this paper, let us list a few physical definitions. One 
possible approach is so-called "internal clock" definition. 
It is based on analysis of LZ probability behavior at finite 
times. As it has been pointed out several times along 
our discussion, the finite time probability dynamics is 
characterized by monotonous function for slow adiabatic 
passage, while for sudden (rapid) transition it oscillates 
before saturation at constant value. These oscillations 
correspond to interference processes and determine the 
population of two states. Therefore, the "internal clock" 
approach defines the Landau-Zener time as the width of 
transition to its asymptotic value (see Refs.fTTl and ITS] 
for detailed discussion). 

An alternative approach to definition the LZ time is 
based on "external clock" probe. In that case the LZ 
Hamiltonian is perturbed by a periodic transverse field 
SH(t) = esin(wt + 4>). Where, to is the frequency of the 
field and cf> its initial phase. The LZ time is determined 
through analysis of infinite time probability as a function 
of external field frequency (see details in Ref . |78j ) . 

Both definitions consistently lead to the estimate of 
LZ times as tlz = A/a for slow adiabatic passage and 
tlz = 1/VO! for rapid passage. Obviously, both defini- 
tions can be straightforwardly generalized for multi-level 
LZ transition. 

Let us consider a slow noise as a special case of "ex- 
ternal clock". We add a perturbation SH = 2f x (t)S x 
to the system such that the coupling A is deviated as 
A{t) — A + f x (t). This case has been discussed in 
the Section [ill B| and corresponds to a non-centered one- 
component transverse noise. The role of noise is to frus- 
trate spin in the direction of the Zeeman field. Let us 
consider a square fluctuation of the Bloch's vector as a 
probe for LZ time: 



{{6b) 2 ) = {b 2 ) - {bf 



(6.1) 



Since the classical noise only dephases the system and 
does not create any dissipation in it, the condition b 2 = 1 
holds. In general, {{8b) 2 ) mixes the diagonal and off- 
diagonal components of the density matrix but for the 
two- and three- level systems subject to classical trans- 
verse noise we consider thus far, {b x ) — {b y ) = and 
b z {t) = pn(t) - f>22{t)- Thus we write, 



{{6b) 2 )=AV sn {t){l-V sn {t)). 



(6.2) 
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FIG. 8. Typical time evolution of the square fluctuations of 
the Bloch's vector. 



The subscript sn refers as usua l to slow noise. It should 
however be noted that relation (6.2 1 works both for two- 
and three- level systems under the assumption that the 



system initially prepared in one of the (upper or lower) 
diabatic states. 

With these ideas in mind, we check the numerical be- 
havior of ((5b) 2 ) for these initial conditions of the spin. 
Essential of results are depicted by Fig[S]below. Interest- 
ingly, ((5b) 2 ) abruptly increases around the anticrossing 
region and saturates to its top value, confirming a spin 
flip transition. After the transition, the variance slightly 
fluctuates (slight decay of ((5b) 2 )) in the direction of the 
Zeeman field for adiabatic addition of noise (see Figs|8b[): 



((5b) 2 ) < ((5b) 2 



(6.3) 



The two- and three- level systems, seems to no longer 
be sensitive to any addition of noise from certain value 
of v. states are in thermal equilibrium. Transition time 
may then be defined as the particular moment where the 
square fluctuation of the Bloch's vector achieved its max- 
imum value. 

This qualitative definition can be experimentally 
probed. For the physical realization of this condition, 
we pose that ((5b) 2 ) is a function of the variable P S n(t)- 
Then, the latter behaves as the quadratic function h(x) = 
4x(l — x) defined in the real space. The mathematical 
requirement for a maximum of a function yields, 



■psn(TLz) 



1 



(6.4) 



Thus, the definition actually coincides with the half- 
width condition V an (rv£) = \V ma x, where V max is the 
maximum value of probability. We emphasize in addition 



that the definition (6.4) holds both for two- and three- 
level systems. 



VII. SUMMARY 

In this paper we discuss the effects of a transverse col- 
ored noise on both two and three-level systems subject to 



LZ transitions. The approximate solution for the tradi- 
tional LZ problem is written down in terms of Fresnel's 
integrals and appears to be useful for exploring the ef- 
fects of slow noise at finite times. We demonstrated in 
the framework of von-Neumann equation, that the ef- 
fects of noise on a two-level system were regulated by an 
integral-differential master equation of the form (2.13|. 



We showed that for a fast Gaussian noise it sufficient 
to average that equation while for the slow noise cor- 
rect procedure is based on averaging the solutions over 
the Gaussian realization of the noise. These arguments 
have been found to be general for description of multi- 
level systems were complicated interference patterns are 
expected. The solution of Bloch's equation is generalized 
for the finite times LZ probabilities of two-and three- 
level system models in the presence of slow noise. We 
have essentially shown that for any initial preparation 
of noise along one of transverse directions (X-noise) or 
two-component transverse noise (XY-noise) the proba- 
bility is renormalized by new functions with shapes of 
standard LZ curves. The famous frequent exponentials 
appearing in LZ transition probabilities are considerably 
discriminated by an inverse square-root function of rf /a. 

In the absence of noise, we showed that population dif- 
ference for two- and three- level system can be found as 
a solution a third order linear differential equation. The 
solution of this equation is given in terms of products 
of the parabolic cylinder Weber's functions. We investi- 
gated solutions by evoking the structural topology which 
generates a logical isomorphism between Schrodinger's 
and Bloch's pictures. In the presence of noise, in general, 
the equation for density matrix elements are integral- 
differential equations. Their solutions can be found 
through the averaging procedure discussed in the paper. 

The appearance of pre-exponential factor depending 
on slow noise in LZ probability for two- and three- level 
systems is consistent with experimental observations re- 
ported in Refs. 44 46]. The suppression of LZ transition 
by the fluctuations of Overhauser field characterized by 
effective temperature associated with noise provides an 
explanation for nuclear spin depolarization mechanism. 



Appendix A: Perturbative solution of Landau-Zener 
problem 



The integral-differential equation for conventional LZ 
problem 



eratively 

5(0) 



2.14) in the absence of noise can be solved it- 



oy setting the perturbative series expansion of 



pW(r) via the parameter A jot: 



P m (r) 



k=0 



(Al) 
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Where, 



x ... 



Pf{r)= f dn r dT 2 CO S [T?-Tl)} 
J — oo J — oo 

/T2k-2 /"T2fe_l 
dr 2 k-i dr 2 feCos[T 2 2 fe _ 1 - T 2 2 fe ], (A2) 

-oo J —oo 



In the presence of noise, we replace A — > 77 in Eq.( |Al[ ). 
Accordingly, the function p^ (r) is modified to p k (r) and 
written as: 



Pfe( T ) = / dT i / dr 2 cos[ri - r|] x ... 

Jo— 2 pT~2h-l 

dT 2 k-l I dT 2k COs[T 2 2 k _ 1 -T 2 2 k \x 



xF«( n ,T 2 ,..,T 2fc ), (A3) 
where the function are defined as follows 

FW(Tl,75,... s 75fc) - 

^ 2fc (/ + (r 1 )/-(T 2 ).../_(T 2fe _ 1 )/ + (r 2fc )). (A4) 

To calculate the higher order correlation function the 
Wick's theorem is used. For the zero-mean random vari- 
ables /+(t) and /_(t) this theorem suggests that: 

F^(r 1 ,r 2 ,...,T 2k ) = V - 2k 

j Ep M r s nLi(/+( T 2«-i)/-(T2n)), for even k, 
1 0, for odd k. 

(A5) 

The summation *^2 pairs runs over all possible combi- 
nations of pairs out of the 2k variables (ri, t 2 , r 2 fc). 
Straightforward calculations for one-component trans- 
verse noise lead to result first derived by Kayanuma^ 

k 

F ik) {T 1 ,T 2 ,...,T 2k ) = ^ eX P (- "f0^2\T 2 „-l ~ T 2n \j, 
pairs n—1 

(A6) 

while for two-components transverse noise (ATK-noise), 
using combinatorics of Keldysh modeP^ we found: 



t 2 , —,r 2 k) 



pairs 



n=l 



^2 2?C ex P ( - 7o ^2\r 2n -i - r 2 

(A7) 



For slow or fast noise, we put respectively 7 — ¥ 0, 7 — > 00 
and the former time-dependent function turns now 
out to be a simple combinatorial factor. The long time 
asymptotic value of P^\t) is evaluated as 



(A8) 



with the help of which we find (2.9| in a perturbative 
expansion form 



a z M = -f>(^)>. 



(A9) 



fc = l 



— l) fe //c! and the radius of convergence for 
Eq. ( A9 1 equal to infinity. Following the argument ac- 



Here, a k 



cording to which for any limit of slow or fast noise the 
probability ( A9 1 is modified by the coefficient in the per- 



turbative expansion, one may write the solution of inte- 
gral equation (3.3) for the cases of slow and fast noise 



driven LZ transition as 



-Plz(oo) = - y^a k b k ( — 
' - \ a 



2 \k 



fe=l 



(A10) 



where b k are the combinatorial factors which appear after 
averaging over realization and depend on either we deal 
with the fast or slow noise. 



1. Slow noise 

It has been shown in Ref.[79] that for slow noise, the 
combinatorial factor b k is expressed as: 



b k =^2 1 = (2fc-l)H, 

pairs 

for an X-noise model and 



b k = J2 2" = 2fefc! > 



(All) 



(A12) 



for an AF-noise. 



pan's 



Fast noise 



The case of fast noise is completely different. In con- 
trast with the slow noise, as shown in Ref.|25| only a 



single term out of the (2k— 1)!! pairings in Eq.(A4) con- 
tributes to the summation for an X- noise: 

b k = 2 k -\ (A13) 

For an XY- model, noise contributes as 

b k = 2 2k -\ (A14) 

3. Discussion 



Equation (A10) can be viewed as result of averaging 



of LZ probability over disorder realizations (the expo- 
nent containing fluctuating field is ensemble average). It 
should be contrasted with the fast noise case for which 
the argument of LZ exponent is proportional to two-point 
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noise correlation function (the argument is disorder aver- 
age). Moreover, the statement concerning disorder aver- 
aging remains true for any finite time transition prob- 
ability. Therefore, the same time dependent function 
F(t) + In W(t) will enter the equation for the finite time 
slow noise driven LZ transition. Since the coefficient bk 
strongly depends on k, the radius of convergence of per- 
turbative expansion must also be found. After summing 
up the perturbative series for LZ probability within the 
circle of convergence, the function has to be analytically 
continued into outer part of the circle. We will identify 
analytical functions describing slow noise driven LZ prob- 



ability and consider finite time LZ transition. Proceed- 
ing, with bk in Eqs.( All )-( A14) we will recover exactly 



all the LZ probability found previously. 



Appendix B: Spin 5 = 1 density matrix equations 

The equation of motion for the density matrix describ- 
ing transitions in there-level system can be obtained the 
same way as we discussed it in details for S = 1/2. In this 
Appendix we present full set of these equations for com- 
pleteness of the discussion about connections between 
Schrodinger and Bloch pictures: 



Here, 



dp+(t) 
dt 

dp-(t) 
dt 



i0+(t)(2pio(t) - A>a(*)) ~ *©-(<) (2poi(<) ~ P2o(t)) 
-»e + (t)(2poa(t) - Pw (t)) + ie_{t)(2p 20 {t) - p 01 (t)). 



Pw(t) = i 



Jl exp (i j£ 6*(r')dT')e_ (*!)/>+ jl e*P (< /£ Q z ^)^ J )& + (t 1 )p li (t 1 )dt 1 



exp 



(t jl®*{T>)dT> 



exp 



(ijl&*(T>)dT> 



Pli{t) = -i 



£ e*P (2^ © Wr') 


e_(ii)(p 2(*i) -ho{h))dh 


exp hi / t * @ z (r')dT^ 





P02(t) = - 



Jl exp (ifl e^dr^Q.^p^dh Jl exp (i j£ 6-(r')dr')e + (t 1 )p 12 (t 1 )di 1 



exp 



(ijlQ*{T>)dT> 



exp 



(iJle*(T>)dT> 



(Bl) 
(B2) 

(B3) 

(B4) 

(B5) 



Here, p %0 {t) = p* t {t) and p±(t) = p* ± (t). 
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